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SECTION  I 


INTRODUCTION  AND  SUMMARY 


1.1  Introduction 

Heat  transfer  from  impinging  jets  is  a  relatively  new  technological 
development  that  has  already  attracted  the  interest  of  heat  transfer  engineers 
and  designers,  involved  with  the  design  of  modern  jet  power  plants  and  related 
machinery.  This  group  of  people  is  mainly  interested  in  heat  transfer  from 
impinging  jets  because  of  its  promise  to  succeed  where  applications  of  convective 
heat  transfer  often  fail  -  in  situations  where  the  rates  of  heat  flux  are 
high  and  the  space  is  restricted  -  such  as  at  the  leading  edge  of  a  turbine 
blade.  This  sudden  spurt  of  interest  of  engineers  primarily  interested  in 
applications  had,  as  a  counterpart,  an  equivalent  response  among  the  more 
research  oriented  circles. 

Investigation  of  the  topics  connected  with  the  heat  transfer  from  the 
impinging  jets  can  be  traced  back  to  about  1951  (Freidman  and  Mueller,  Reference 
1)  for  an  experimental  approach  and  to  1956  (Kezios,  Reference  2)  for  a 
combination  of  analysis  and  experimentation.  Since  that  time,  many  investigators 
have  appeared  who  tried  both  analytical  and  practical  approaches,  mostly 
with  one  ultimate  purpose  in  mind:  to  furnish  to  the  designer  formulas,  charts 
and  graphs,  that  would  enable  him  to  handle  various  heat  transfer  problems 
related  to  the  impinging  jets  in  a  routine  fashion. 

In  general,  the  geometry  of  even  a  single  impinging  jet  is  quite  complicated, 
and  the  heat  transfer  aspects  depend  strongly  on  the  corresponding  fluid 
flow  structure.  Therefore,  in  the  following  literature  review,  fluid  flow 
aspects  of  the  impinging  jets  will  be  reviewed  first,  and  then  our  attention 
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will  be  turned  to  the  items  more  directly  related  to  heat  transfer  from  the 
Impinging  jets. 

Among  the  other  hydrodynamic  parameters,  one  of  the  most  critical  parameters 
governing  heat  transfer  from  Impinging  jets  appears  to  be  the  distance  between 
the  nozzle  exit  and  the  target.  In  particular,  the  heat  transfer  from  the 
impinging  jets  to  flat  targets,  appears  at  first  to  Increase  directly  with  the 
distance  between  the  nozzle  and  the  target,  to  reach  a  flat  maximum  at  six  or 
seven  diameters  away  in  the  case  of  circular  jets.  This  distance  Is  for  the 
slot  jets  slightly  higher,  and  for  both  kinds  of  jets,  the  region  of  maximum 
intensity  in  heat  transfer  coincides  roughly  with  the  tip  of  the  "potential 
core"  -  the  region  where  the  original  jet  velocity  at  the  nozzle  exit  is  still 
preserved  at  the  center  of  the  jet. 

Beyond  the  potential  core  region,  the  intensity  of  heat  transfer  decreases 
pari  passu  with  the  further  dissipation  of  the  original  jet  velocity.  This 
zone  of  jet  development  has  originally  attracted  a  substantial  number  of  in¬ 
vestigators  who  were  concerned  primarily  with  the  free  jets,  unencumbered  by 
the  presence  of  the  targets.  Here,  similarity  profiles  techniques  have  been 
used  with  success. 

After  the  jet  hits  the  target,  eventually  a  wall  jet  develops,  which  ex- 
nlbits  definite  boundary  layer  characteristics.  Consequently,  within  the 
potential  core,  heat  transfer  intensity  will  depend  on  the  original  jet  velocity, 
whereas  in  the  region  of  the  developing  jets,  heat  transfer  Intensity  is  a 
function  of  the  actual  velocity  of  the  jet  when  it  hits  the  target  the  so 
called  "velocity  of  arrival". 
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1.2  Summary 

In  the  report,  the  available  literature  on  the  subject  of  flow  patterns 
and  heat  transfer  resulting  from  impinging  jets  has  been  reviewed.  The 
report  shows  that  despite  some  accomplishments,  there  still  exists  necessity 
for  both  analytical  and  experimental  work,  before  the  heat  transfer  aspects 
of  the  impinging  jets  are  properly  understood.  The  state  of  the  art  with 
respect  to  fluid  flow  patterns  is  in  a  relatively  better  position,  but  there, 
too,  appear  considerable  gaps  in  knowledge,  primarily  concerning  the  effects 
of  turbulence  intensity  and  due  to  a  diff icult-to-handle  geometry. 

The  flow  from  free  jets  becomes  turbulent  in  the  Reynolds  number  range 
from  1000  to  3000.  Impinging  jets  behave  still  like  free  jets  with  the 
target  only  slightly  over  one  nozzle  diameter  away  from  the  jet  nozzle. 

The  nominal  potential  core  region,  for  the  non-turbulized  jets  only, 
extends  six  to  seven  diameters  away  from  the  nozzle  for  the  circular  jets; 
for  slot  jets,  this  distance  seems  to  be  somewhat  longer. 

The  flow  near  the  stagnation  point  is  laminar,  but  the  effects  of  the 
jet  stream  turbulence  on  the  heat  transfer  after  impingement  are  considerably 
higher  than  for  the  flow  over  a  flat  plate.  This  makes  it  sometimes  difficult 
to  determine  how  far  away  from  the  stagnation  point  the  laminar  regime 
still  extends.  Outside  boundary  layer,  potential  flow  solutions  are  valid 
near  the  stagnation  point. 

Because  in  the  real  jets,  of  a  finite  thickness,  the  flow  regime  depends 
considerably  on  the  distance  away  from  the  nozzle,  it  seems,  at  the  present 
time,  that  with  a  few  exceptions,  analytical  calculations  are  only  possible 
when  the  target  is  relatively  close  to  its  nozzle  so  that  flow  may  still  be 
considered  infinite,  in  the  sense  as  originally  treated  by  Hiemenz  and  Homann. 

However,  there  appears  to  be  a  real  promise  in  semi-empirical  calculations 
of  the  hydrodynamic  and  heat  transfer  phenomena  in  the  wall  jet  region, 
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farther  away  from  the  stagnation  point,  where  the  flow  Is  alreay  definitively 
turbulent.  For  arrays  of  jets,  the  results  so  far  have  been  obtained  only 
experimentally,  and  there  still  exist  gaps  in  our  knowledge  of  the  subject. 

It  seems  that,  at  the  present  state  of  the  art,  the  calculations  of  heat 
transfer  from  arrays  of  jets  will  have  to  rely  heavily  on  experimental  techniques, 
or  on  analogue  computation  for  some  time  to  come. 

It  Is  interesting  to  note  the  effects  of  the  turbulence  intensity  and 
the  turbulence  scale  on  local  heat  transfer  near  the  stagnation  point,  for 
both  the  two-and  the  three-dimensional  jets,  impinging  on  flat  targets. 

These  effects,  observed  independently  from  each  other  by  several  investigators, 
manifest  themselves  in  a  substantial  increase  of  heat  transfer  intensity 
in  the  neighborhood  of  the  stagnation  point,  and  in  an  explicit  dependence 
of  heat  transfer  intensity  there  on  the  no2zle  diameter. 
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SECTION  II 

FLUID  FLOW  ASPECTS  OF  IMPINGING  JETS 
2.1  General  Features  of  Jets,  Potential  Core  Length 

As  in  virtually  all  forced  convection  problems,  fluid  flow  aspects  of 
impinging  jets  cannot  be  really  separated  from  the  heat  transfer  aspects 
(and  vice  versa),  since  the  energy  equation  is  dependent  upon  the  momentum 
equation.  A  reasonable  degree  of  understanding  of  fluid  flow  patterns  associated 
with  the  impinging  jets,  is  the  prerequisite  for  a  satisfactory  analysis  of 
all  the  associated  heat  transfer  effects. 

Of  fundamental  importance  is  here  the  distinction  between  the  laminar 
and  the  turbulent  flow  regimes.  Here,  according  to  Vickers,  (Reference  3), 
it  seems  that  the  critical  Reynolds  number.  Re,  based  on  the  nozzle  diameter, 
is  most  likely  1000.  This  finding  is  supported  by  McNaughton  and  Sinclair 
(Reference  4),  who  report  having  observed  four  characteristic  jet  patterns: 

1.  What  they  call  a  "dissipated  laminar  jet",  Re  <  300. 

2.  Fully  laminar  jet,  300  <  Re  <  1000. 

3.  Semi -turbulent  jet,  that  starts  as  a  laminar  jet,  eventually 
becoming  turbulent,  1000  <  Re  <  3000. 

4.  Fu'ly  turbulent  jet.  Re  >  3000. 

The  above  observations  were  made  on  liquid  into  liquid  jets,  in  vessels 
of  various  shapes,  using  as  a  tracer  methylene-blue  dye.  Working  with  air 
jets,  Hrycak  et  al.  (Reference  5)  have  obtained  very  similar  results. 

The  above  results  are  in  line  with  the  observations  made  by  Cederwall, 
(Reference  6)  who  found  that  jets  with  Re  >  3000  are  already  turbulent.  More¬ 
over  the  heat  transfer  results  of  Gardon  and  Cobonpue,  (Reference  7)  and 
those  by  Chamberlain,  (Reference  8)  seem  to  indicate  that,  for  impingement  in 
the  potential  core  region,  the  Nusselt  number  ceases  to  be,  (for  the 
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given  fluid  and  flow  geometry)  a  function  of  the  Reynolds  number  alone  for 
Re  >  7000.  This  means  that  the  Nusselt  number  must  become  then  also  functionally 
dependent  on  turbulence  patterns  prevailing  in  the  jet,  as  the  simplest 
explanation  of  this  phenomenon. 

The  above  remarks  apply  directly  only  to  jets  issuing  from  circular 
openings,  but  are  also  roughly  applicable  to  slot  jets.  Thus,  in,  another 
paper  dealing  with  slot  jets  exclusively,  Gardon  and  Akfirat,  (Reference  9) 
consider  jets  above  Re  =  2000  as  turbulent. 

The  above  discussion  applies,  strictly  speaking,  only  to  free  jets.  It 
seems,  however,  that  similar  reasoning  can  be  applied  also  reasonably  well 
to  the  impinging  jets,  in  the  space  not  affected  by  the  presence  of  the 
target  plate.  Theoretical  calculations  show  that  the  effect  of  the  impinging 
surface  is  felt  by  the  jet  only  in  the  immediate  vicinity  of  that  surface 
(Levey,  Reference  10).  Even  when  the  impinging  surface  is  only  one  slot 
width  away  from  the  origin  of  the  jet,  it  has  no  effect  on  either  the  momentum 
flux  or  the  mass  flux  through  the  jet  in  question.  For  axi symmetric  jets 
this  result  is  essentially  confirmed  by  Tani  and  Komatsu’s  experimental 
observations,  (Reference  11)  and  Hrycak  et  al.  (Reference  5).  It  should  be 
noted  that  the  above  observed  critical  Reynolds  number  values  are  one  order 
of  magnitude  higher  than  the  results  obtained  from  the  flow  stability  calcula¬ 
tions,  as  outlined  in  Reference  12,  for  example. 

Another  important  characteristic  of  a  jet  is  the  length  of  the  potential 
core  (Figure  1)  -  the  distance  to  which  the  original  nozzle  exit  velocity  of 
the  jet  persists  along  the  center  line  of  the  jet. 

Here,  some  useful  theory  is  available.  According  to  Schlichting,  (Reference 
12,  p.  732),  we  have  the  width  of  the  jet,  b,  proportional  to  the  distance 
from  the  nozzle  exit,  x. 

b  =  const,  x  (1) 
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The  relationship  between  the  maximum  velocity,  u  ,  at  the  jet  centerline, 
and  x,  can  be  obtained  from  the  conservation  of  momentum  equation,  where  the 
x-component  of  total  jet  momentum,  J,  must  be  independent  of  x: 

J  =  /  pu2  dA  =  const.  (2) 

A 

2 

In  the  case  of  two-dimensional  jet,  we  have  J*  =  const,  x  pu^b  where  J'  * 
momentum  per  unit  length.  Equations  (1)  and  (2)  may  be  used  to  get  expressions 
for  um>  valid  in  the  region  where  um/uoc  <  1  applies.  In  the  equations  for  um 
below,  the  constants  and  C2  must  be  determined  empirically  for  the  best  re¬ 
sults.  Stated  In  dimensionless  form,  we  have  then  for  um/uQC  the  relations 


un/uoc  =  C1  '  Hl\  *  xo/Bo 


for  the  slot  jet,  and 


(3) 


V“oc  *  C2/(x/D  *  v°>  (4) 

for  the  circular  jet.  Equations  (3)  and  (4)  are  valid  only  outside  of  the 
potential  core  region;  xQ  term  Is  the  virtual  origin  correction,  since  a  real 
Jet  does  not  start  from  a  point  source,  but  at  the  nozzle  exit  with  a  finite 
diameter  and  finite  velocity  uQC  at  the  center  of  jet. 

The  above  results  can  be  used  to  determine  the  length  of  the  potential 
core  region,  when  the  ratio  u  /u  is  plotted  against  the  dimensionless 
distance  from  the  jet  nozzle  exit  on  log-log  paper.  The  slope  of  the  resulting 
curves  Is  Indeed  very  nearly  -1/2  for  the  slot  jets,  and  -1  for  the  circular 
jets  (Reference  5).  The  values  of  the  constants  and  C2  determine  the 
length  of  the  nominal  potential  core  distance.  The  actual  values  of  these 
constants  depend  on  the  skill  of  the  particular  investigator,  the  specific 
method  used  to  get  and  C2,  as  well  on  the  nozzle  construction  details. 

For  example,  the  value  of  2.28  for  C1  has  been  recommended  by  Albertson  et 
al.  (Reference  13)  which  makes  the  end  of  the  potential  core  5.2  slot  widths 
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away  from  the  slot  Itself,  whereas  Schauer  and  Eustis  (Reference  14)  used 
the  value  of  C i =2 . 35  according  to  the  results  observed  by  Van  der  Hegge  Zljnen 
(Reference  15). 

For  Cj,  Albertson  et  al.,  give  the  value  of  6.2.  In  discussion  on 
Reference  13,  Citrini  recommends  an  empirical  relationship 

u/u  =  6.6  D/x  -  0.49 

m  o 

while  Baines  claims  that  there  Is  a  lengthening  of  the  potential  core  with 
Increasing  Reynolds  number  from  5.0  at  Re  s  1.4x10^,  to  7.0  at  Re  *  1.0x10^ 
Abramovich  (Reference  16),  on  the  other  hand,  asserts  that.  C2  does  not  vary 
with  the  Reynolds  number  directly,  but  does  so  with  the  velocity  profile. 

Since  the  velocity  profile  in  turbulent  flow  is  a  weak  function  of  the 
Reynolds  number,  it  means  that  C2  varies  (although  to  a  very  limited  extent) 
also  with  the  Reynolds  number.  After  all,  the  recorded  variation  of  C2  Is 
from  5.2  (adopted  from  Tani  and  Komatsu  Reference  11)  to  7.7  (Poreh  and 
Germak,  Reference  17);  all  these  figures  are  based  on  the  average  velocity  at 

the  jet  nozzle  exit.  The  mean  value  of  C2  may  be  taken  as  7  (cf.  Gauntner  et 
al..  Reference  18).  A  comprehensive  study  of  the  magnitude  of  C2  has  been 
carried  out  by  Hrycak  et  al.  (Reference  5). 

Another  aspect  of  the  potential  core  length  is  its  variation  with  the 
Mach  number  as  discussed,  for  example,  by  Snedeker  and  Donaldson  (Reference 
19)  who  find  an  apparent  potential  core  length  of  7.5  diameters  for  M  = 

0.52.  Downstream  of x  /D  *11,  the  velocity  decay  follows  the  1/x  dependence 
quite  closely.  Here  Snedeker  and  Donaldson  cite  the  experiments  of  Warren 
(Reference  20),  who  for  M  «  0.69  found  a  core  length  of  7.2  diameters  with 
x/D  c  10  for  the  start  of  the  fully  developed  region. 

Actually,  the  exact  "measured"  length  of  the  potential  core  depends  on 
the  method  of  Its  determination.  Since  there  is  no  sudden,  sharp  transition 
at  the  center  between  the  potential  core  and  the  region  of  the  developed 


(5) 
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flow,  this  determination  must  be,  by  necessity,  a  somewhat  arbitrary  procedure. 
Thus,  for  example,  Pal  (Reference  21)  suggests  for  the  constant  in  Equation 
(4)  the  value  =  6.5. 

Pai  reports  also  the  existence  of  the  intermittancy  in  the  jet  flow, 
first  observed  by  Corrsin:  in  a  fully  developed  axi -symmetric  jet,  a  completely 
turbulent  flow  exists  only  in  the  core  region,  out  to  a  radius  at  which  the 
velocity  is  about  one  half  of  the  maximum  velocity  of  the  given  cross- 
section.  Outside  of  this  region,  there  is  a  transition  region,  followed  by 
a  region  of  laminar  flow. 

Information  about  the  exact  length  of  the  potential  core  is  quite  important 
in  regard  to  maximum  cooling  efficiency;  heat  transfer  results,  obtained  with 
single  jets,  show  that  the  maximum  intensity  of  heat  transfer  takes  place  at 
the  tip  of  the  nominal  potential  core. 

The  above  results  on  the  form  of  attenuation  of  the  jet  centerline 
velocity  with  the  distance  can  also  be  calculated  through  assumption  of  a 
particular  similar  velocity  profile,  with  a  specific  velocity  distribution, 
for  flow  regions  where  similar  velocity  profiles  can  be  expected  to  have 
already  established  themselves  (Albertson  et  al.,  Reference  13),  and  then 
working  backwards,  towards  u/u  =  1.0..  (This  is  of  course  only  possible 

vV 

when  the  potential  core  has  already  disappeared,  and  should  be  applicable, 

In  principle,  to  both  the  laminar  and  the  turbulent  flow  regimes). 

There  exist  essentially  two  methods  of  determining  the  potential  core 
length  experimentally.  The  real  potential  core,  L*,  ends  where  the  jet 
centerline  velocity  starts  to  deviate  for  the  first  time  from  uQC.  The 
"nominal"  potential  core  length  is  obtained  from  an  extrapolation  of  the 
developed  jet  centerline  velocities  plot  on  log-log  paper  towards  lower  x/B0 
or  x/D  values.  Equation  (3)  and  (4)  here  refer  to  the  nominal  potential 


core  length;  the  nominal  potential  core  length  is  a  useful  concept  in  heat 
transfer  calculations,  however.  It  is  about  2.3  nozzle  diameters  longer 
than  the  true  potential  core,  L*. 

In  general,  in  the  case  of  an  impinging  jet  we  can  distinguish  four 
characteristic  regions  of  flow  (Figure  1).  They  are  (Gauntner  et  al..  Reference 
18): 

1.  a  transition  zone  of  flow  establishment 

2.  a  zone  of  established  flow  in  the  original  direction  of  the  jet 

3.  a  deflection  zone 

4.  a  zone  of  established  flow  in  the  radial  direction. 

In  Zone  1,  the  most  interesting  fact  is  the  existence  of  the  potential 
core  region  with  the  undisturbed  original  velocity  of  the  jet,  already  discussed 
above.  Ordinarily,  Zone  1  extends  up  to  nine  nozzle  diameters  downstream 
from  the  jet  orifice.  In  Zone  2,  the  velocity  distribution  may  be  considered 
similar  to  that  of  a  free  jet  diffusing  into  an  infinite  medium.  In  the 
neighborhood  of  the  stagnation  point,  (Zone  3),  flow  patterns  appear  to  be 
close  to  those  computed  for  the  irrotational  flow. 

Of  great  importance  is  the  zone  of  established  flow  in  the  radial  direction. 
Zone  4.  The  jet  in  this  region  is  called  the  wall  jet,  and,  for  axisymmetric 
jets,  exhibits  characteristics  somewhat  similar  to  those  of  a  free  radial 
jet  -  a  jet  issued  from  a  toroidal  slit  in  the  radial  direction.  The  presence 
of  the  impinging  plate  is  responsible  for  the  boundary-layer  development. 

In  many  practical  applications,  Zone  2  of  jet  flow  has  never  a  chance 
to  form  if  Zn/D  <  9,  Zn  being  the  normal  distance  between  the  target  plate 
and  the  nozzle.  Zone  \  has  been  already  discussed  to  some  extent.  In  the 
following,  the  Implications  of  existence  of  Zones  3  and  4  will  be  covered. 

In  Zones  3  and  4.  th"  boundary- layer  theory  can  be  applied  with  some 
far-rrarhing  r.nnsrfjurnces .  |t  is  particularly  convenient,  that  near  the 
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stagnation  point,  the  available  exact  analytical  solutions  of  the  complete 
Navler  -  Stokes  equations,  applicable  for  the  laminar  flow  Impinging  on  a  wedge, 
can  be  utilized  (Reference  12).  Brady  and  Ludwig  (Reference  22)  prepared 
a  chart  to  Indicate  the  extent  of  the  laminar  region  near  the  stagnation  point. 

A  summary  of  fluid  flow  characteristics  of  a  round  turbulent  jet  Impinging  on 
a  flat  plate  has  been  made  by  Gauntner  et  al..  Reference  18. 

2.2  Analytical  Treatment  of  Impinging  Jets 

The  simplest  approach  to  the  analytical  handling  of  impinging  jets  is  based 
around  the  assumption  of  inviscld  flow.  The  essential  difficulty  here  lies  in 
the  fact  that  despite  the  well  defined  boundary  conditions  in  the  form  of  the 
free  streamlines  in  such  flows,  the  location  of  the  free  streamlines  Is  not 
known  in  advance,  and  must  be  determined  as  a  part  of  the  solution  (Reference  23). 

Exact  analytical  solutions  for  two-dimensional  Impinging  flow  are  available 
for  Ideal  fluids  (e.g.  Levey,  Reference  10).  The  available  theoretical  treatments 
of  Inviscld  three-dimensional  Impinging  jets  have  been  obtained  through  analog 
solutions  or  approximate  analytical  and  numerical  methods  (e.g..  Strand,  Schach, 
LeClerc,  and  Dal bes, References  24-27).  The  results  of  the  researchers  above  are 
sound,  because  they  check  quite  well  with  the  available  experimental  data.  The 
approach  based  on  the  Inviscld  flow  can  be  expected  to  give  solutions  which  are 
reasonably  close  to  an  actual  situation  where  the  dynamic  effects  predominate  over 
the  viscous  effects,  which  occurs  at  sufficiently  high  Reynolds  numbers.  Here, 
Schach  used  an  integral  equation  method  according  to  Trefftz,  which  resulted  In 
numerical  solutions  by  successive  approximations,  for  the  location  of  the  free 
jet  boundary.  LeClerc  used  an  electrolytic  analog  to  locate  the  free  jet 
boundary  and  to  determine  the  shape  of  the  streamlines. 


In  his  approach,  once  the  free  surface  of  the  jet  had  been  obtained, 

LeClerc  employed  relaxation  techniques  in  the  interior  of  the  jet  to  establish 
the  potential  lines  and  streamlines.  Daibes  obtained  an  approximate  solution 
for  an  Invlscld,  Impinging  jet  using  cylindrical  coordinates  to  determine  the 
free  surface  shape  from  consideration  of  the  exact  boundary  conditions.  The 
governing  differential  equations  were  then  solved  by  relaxation  methods. 

Daibes'  numerical  approach  essentially  confirmed  Schach’s  and  LeClerc's 
analyses  (Figure  2). 

The  velocity  distribution  (Figure  3)  Is  connected  by  Equation  (6) 

v/uQ  =  /  1  -  (p-pQ)/*5  p  u*  (6) 

to  the  pressure  distribution.  At  the  stagnation  point,  potential  flow  solutions 
result  In  v/u  =  a*  r/D,  with  a*  =  1.02  by  Daibes,  and  a*  =  1.03  by  LeClerc. 

The  dimensionless  velocity  gradient,  a*.  Is  needed  for  heat  transfer  calculations 
at  the  stagnation  point,  and  will  be  the  subject  of  additional  discussion  below. 

In  the  region  where  the  laminar  flow  prevails,  exact  solutions  of  the 
Navler  -  Stokes  equations  are  possible  (Reference  12).  It  also  happens  to  be 
the  case  here  that  the  boundary  layer  approach  gives  the  same  results  as  the 
corresponding  solutions  of  the  complete  Navler  -  Stokes  equations.  The  flow 
within  the  boundary  is  dependent  upon  the  main  flow,  in  the  sense  that  the 
flow  configuration  determines  the  pressure  distribution  at  the  surface  under 
consideration  (References  2,  12). 

For  systems  of  axial  symmetry,  momentum,  energy  and  continuity  equations 
are  then, 

p(v  If  +  u  +  ll  If)  (amentum)  (7) 

where  for  z-0:  u=0,  v=0;  for  z=»  :  v°vs 

cP»t#!?*U!l>  B5F(kll>  <8> 
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where  for  z=0,  T=TW;  for  z=“  ,  T=T» 

I?  rv^  +  fjr  (p  ru)  =  0  (continuity)  (9) 

When  the  free  stream  velocity  is  proportional  to  some  power  of  the  dis¬ 
tance  from  the  stagnation  point,  this  leads  to  the  notion  that  the  heat  transfer 
coefficient  at  any  point  along  the  median  profile  of  a  body  of  arbitrary  cross- 
section  is  the  same  as  that  on  a  wedge,  which  has  the  same  velocity  and  pressure 
gradient  just  outside  the  boundary  layer  at  the  same  distance  from  the  stagnation 
point  as  the  real  profile.  In  addition,  for  heat  transfer  applications  the  local 
temperature  ratio  has  to  be  the  same  (Reference  2). 

Thus,  according  to  this  procedure,  the  two-dimensional  surface  is  replaced 
by  a  series  of  wedges,  each  having  a  uniquely  enclosed  angle  that  depends  on 
the  Euler  number,  m.  In  question.  The  three-dimensional  stagnation  flow  repre¬ 
sents  a  special  case  of  wedge  flow. 

In  the  stagnation  region  (Zone  3  in  Figure  1)  the  velocity  at  stagnation 
point  obeys  the  relation  v  *  ar;  when  a*  =  aD/u_r  is  Introduced,  one  has  outside 
of  the  boundary  layer  also  the  relation  v  =  u  a*r/D;  the  system  of  Equations 

v  V/V* 

(7)  -  (9)  can  be  solved  (Reference  12)  by  reducing  it  to  the  total  differential 

equation  in  f(n),  by  using  n  *  v^Tvz.t  =  r2  v'ivf,  vr  =  and  ur  =  -  —  . 

qZ  of* 

Letting  f'  *  df/dn,  etc.,  there  results 

f"  +  2  ff"  +  1  -  f’2  *  0  (10) 

where  f  =  0,  f'  =  0  at  n=  0,  f*  ■  latn=  »,  and  v/v$  «  f'.  This  kind  of  solution 

was  first  obtained  by  Homann  (see  Reference  12,  p.  100).  Similarly,  one  can 

define  a  dimensionless  temperature  function,  (T-T  )/(T  -T  )  =  g(n);  through 
substituting  for  T,  u  and  v  In  terms  of  f  and  g  and  their  derivatives,  one 
gets  from  the  energy  equation  the  dimensionless  expression 

g"  +  2  f  Pr  g'  =0 
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where  g  s  1  atn*  0,  and  g  ■  0  at  n  s  •.  Solution  of  Equation  (11)  has  been 
obtained  by  Sibulkln  (Reference  28) 

Nup/Re^  -  0.763  Pr0,4  (a*)1*  (12) 

where  a*  Is  a  free  constant  that  depends  on  a  particular  system's  geometry: 
for  example,  a*  *  3  for  a  sphere,  and  a*  -*  1.0  as  Zn/D  -*•  1 ,  In  the  limiting 
cases  of  impinging  jet  flow. 

From  looking  on  Equations  (9)-(12)  it  is  seen  that  effects  of  fluid  flow 
are  intimately  interrelated  with  the  heat  transfer  aspects  of  the  impinging 
jets.  Thus,  the  Nusselt  number  from  Equation  (12)  is  a  function  of  the 
Reynolds  number  that  is  linked  to  the  flow  regime  and  the  flow  geometry,  and 
to  the  Prandtl  number,  which  brings  in  the  thermal  properties  that  determine  the 
physical  aspects  of  heat  transfer.  Finally,  there  Is  the  a*  term,  the  dimension¬ 
less  velocity  gradient  at  the  stagnation  point;  a*  is  the  dimensionless  ratio 
that  relates  the  fluid  flow  regime  at  the  stagnation  point  to  the  conditions  at 
the  nozzle  exit  and  to  the  nozzle-to-target  distance. 

In  practical  cases,  a*  is  a  function  of  the  jet  nozzle  velocity  profile  at 
the  exit  (which  Is,  In  general,  also  related  to  the  nozzle  geometry)  and  the 
distance  between  the  target  plate  and  the  nozzle  exit,  In  the  case  of  real  fluids. 
Turbulent  mixing  at  the  boundary  of  the  potential  core  (for  ReD  >  1000)  modifies 
the  original  velocity  profile  as  a  function  of  the  distance  from  the  nozzle;  con¬ 
sequently,  a*-magnitude  generally  decreases  with  the  distance  between  the  nozzle 
and  the  target  plate,  In  real  life  situations  (cf.  Reference  5).  Additional 
discussion  of  a*  will  be  found  In  Section  IV,  dealing  with  experimental  Impinging 
jet  heat  transfer. 

Stagnation  flow  changes  eventually  Into  wall -jet  flow  (Zone  4  in  Figure  1) 
after  the  direction  of  flow  has  been  turned  through  a  90°  angle.  In  the  wall- 
jet  flow  of  real  fluids,  we  may  distinguish  the  outer  region,  z  >  6,  where  the 
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effects  of  friction  between  the  fluid  in  the  jet  and  the  quiescent  surrounding 
fluid  are  predominant,  and  the  Inner  region,  z<6,  where  the  presence  of  the  target 
plate  makes  itself  felt.  Because  each  region  has  distinct  physical  character¬ 
istics,  the  flow  is  only  approximately  similar  in  Zone  4. 

The  flow  In  Zone  4  has  been  analyzed  by  Glauert  (Reference  29).  Although 
some  of  the  assumptions  used  in  that  analysis  proved  to  be  only  remotely  related 
to  reality  (for  example,  the  existence  of  wall  friction  coefficient  of  the  Blasius 
type),  Glauert's  solution  describes  reasonably  well  flow  situation  in  the  wall- 
jet  zone,  and  is  potentially  useful  for  the  analysis  of  heat  transfer  there, 
with  certain  limitations.  Also,  heat  transfer  in  Zone  4  may  be  still  analyzed 
by  semi -empirical  methods.  One  such  method  is  based  on  the  Colburn  analogy, 
applicable  to  fluids  whose  Prandtl  number  is  sufficiently  close  to  unity: 

St  Pr2/3  =  Cf/2  (13) 

where  the  Stanton  number  is  defined  as  St  =  Nu/Re  Pr  =  h/(vCp).  If  the  friction 

coefficient  Cf  is  known,  as  a  function  of  the  local  Reynolds  number,  and  ether 

pertinent  wall  jet  parameters  are  available,  the  heat  transfer  coefficients 
(local  and  average)  can  then  be  calculated.  For  a  turbulent  wall  jet,  the  use 
of  Equation  (13)  appears  to  be  the  most  convenient  way  of  calculating  the  wall- 
jet  heat  transfer. 

Stagnation  point  mass  transfer  follows  governing  equations  similar  to  those 
for  heat  transfer.  For  identical  boundary  conditions,  Nu/(Pr)n  =  Sh/(Sc)n, 
where  Sh  =  Sherwood  number,  and  Sc  =  Schmidt  number,  are  the  Nusselt  number, 
and  Prandtl  number,  respectively,  equivalents  for  mass  transfer. 

Many  other  deserving  authors,  who  have  written  on  the  subject  of  fluid 
flow  patterns  of  free  or  impinging  jets,  have  not  been  mentioned  here  directly. 
This  was  done  so  because  of  the  time  constraints  Imposed  upon  the  present 
study,  supposed  to  serve  mainly  as  an  Introduction  to  the  bibliography,  rather 
than  as  an  exhaustive  survey  article,  concerning  the  impinging  jets. 
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However,  It  Is  hoped  that  the  discussion  above  has  elucidated  some  of 
the  facts  that  make  analytical  calculations  of  fluid  flow  parameters  and,  con 
sequently,  of  heat  transfer  from  Impinging  jets,  so  difficult.  This  Is  so 
because  of  considerably  more  Involved  geometry,  and  more  highly  pronounced 
frictional  effects.  In  comparison  with  the  corresponding  analytical  and  seml- 
emplrlcal  solutions  available  for  a  flat  plate  In  the  parallel  flow,  or  In 
the  generalized  stagnation  flow  (not  restricted  to  a  relatively  thin  jet). 


SECTION  III 

ANALYTICAL  HEAT  TRANSFER 

3.1  Results  Obtained  with  Standard  Analytical  Approaches 

The  region  of  the  jet  impingement.  Zone  3,  of  Figure  1  was  investigated 
by  numerous  individuals  who  employed  similarity  techniques  adapted  to  stagnation 
flow  and  the  wedge  -  flow  solutions  of  the  boundary  layer  equations.  After 
the  jet  hits  the  target,  it  is  deflected  and  eventually  a  wall  jet  develops 
(Zone  4,  of  Figure  1).  Since  the  wall-jet  region  exhibits  boundary- layer 
characteristics,  it  may  be  analyzed  by  boundary-layer  theory. 

Initial  attempts  to  analytically  predict  heat-transfer  from  an  impinging 
jet  to  a  plane  surface  were  limited  to  the  neighborhood  of  the  stagnation 
point,  where  even  for  turbulent  jets,  laminar  flow  predominates  in  the  boundary 
layer.  Kezios,  (Reference  2)  reports  the  following  expression  for  the 
stagnation-point  heat  transfer: 

Nuo  =  0.475  /KS^  (14) 

where  the  nozzle  radius  D/2  was  used  as  the  characteristic  dimension.  Equation 
(14)  becomes,  when  Reg  is  used  and  the  term  Pr0>4  is  introduced,  identical 
with  Sibulkin's  solution.  Kezios  found  also  analytically,  at  r/D  -  0.65,  a 
maximum  for  the  local  Nusselt  number,  away  from  the  stagnation  point. 

Walz  (Reference  30)  using  the  results  originally  derived  by  Eckert,  obtained 
the  expression 

Nu  =  0.44  Pr°-36Rer0’5  (15) 
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for  the  local  values  in  the  axi symmetrical  flow,  with  the  average  value  of 
the  Nusselt  number  in  the  axi symmetrical  flow  being  equal  to  twice  the  local 
value. 

For  slot  jets,  using  Eckert's  wedge  flow  solutions,  Metzger  (Reference 
31)  obtained  the  equation 

Nu  *  0.57  Pr0,37  Rer0,5  (16) 

for  local  values  of  the  Nusselt  number  based  on  distance  from  the  stagnation 
point,  and  found  the  average  value  of  the  Nusselt  number  to  be  twice  the 
local  value.  The  inability  of  such  formulas  to  show  explicitely  the  effect  of 
the  nozzle-to-plate  separation  distance  may  be  considered  a  weakness. 

Schauer  and  Eustis  (Reference  14)  present  the  solution  for  the  hydrodynamics 
of  a  two-dimensional  wall  jet,  the  basis  for  the  analysis  in  Reference  32. 

In  Reference  32,  the  thermal  boundary- layer  equations  for  incompressible 
turbulent  flow  were  solved  by  applying  integral  techniques.  A  modified 
Reynolds  analogy  (relating  the  diffusivity  for  heat  to  the  diffusivity  for 
momentum)  was  used  to  link  the  thermodynamic  to  the  hydrodynamic  solution. 

In  the  impingement  region,  of  a  two-dimensional  jet,  Cadek  (Reference 
33)  used  an  approach  that,  although  based  on  the  integral  techniques,  was 
reducible  to  Eckert's  wedge  flow  solutions.  In  the  wall  jet  region,  Cadek’s 
efforts  resulted  in  expressions  which  could  be  reduced  to  the  well-known 
flat-plate  solution  for  the  corresponding  limiting  case. 

Cadek  recognized  the  fact  that, in  the  wedge  flow  class  of  solutions,  an 
infinite  uniform  stream  impinges  on  a  surface,  whereas  in  the  case  of  impinging 
jets  we  deal  with  a  nonuniform  jet  of  finite  width  striking  a  surface.  In 
the  first  case,  the  Euler  number  m  is  constant;  in  the  latter  case,  it  is  a 
variable  function  of  distance  from  the  stagnation  point. 
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Eckert  and  livingood  developed  methods  of  calculating  heat  transfer  with  a 
variable  m  in  Reference  34,  applicable  to  stagnation  flow.  According  to 
Cadek,  his  method  gives  results  on  heat-transfer  rates  that  were  essentially 
similar  to  those  obtainable  by  Eckert  and  Livingood' s  approach,  while  at 
the  same  time  yielding  more  information  on  the  boundary  layer  itself. 

Tomich  calculated  heat  transfer  from  impinging,  turbulent  jets  (Reference 
35)  by  solving  the  governing  equations  by  means  of  finite  difference  techniques 
that  were  extended  into  the  compressible  flow  region.  !t  was  found  that  the 
jet  Mach  number,  and  jet  temperature  ratio,  were  the  only  two  initial  jet 
properties  necessary  to  characterize  the  dimensionless  velocity  and  temperature 
variations  in  the  free  jet. 

In  Reference  36,  Brdlik  and  Savin  used  the  integrated  energy  equation 
adapted  for  situations  with  a  radial  symmetry  and  obtained,  for  Zn/D<6.0,  with  Re 
based  on  the  plate  radius  Dp/2 

No  =  1.09  Pr1/3Re0p/2  U7) 

which  shows  reasonably  good  agreement  with  the  test  data  cf  several  invest  gators 
for  0.3  <  Pr  <  1.7. 

Sparrow  and  Wong  derived,  (Reference  37)  analytical  expressions  for 
heat  transfer  from  a  slot  jet  with  a  parabolic  velocity  profile  that  shows  a 
26%  increase  for  Zn/B0=l,a  68%  increase  for  zn/B0=2,  etc.,  in  comparison 
with  heat  transfer  from  a  slot  jet  with  a  flat  initial  velocity  profile  at 
the  nozzle. 

Miyazaki  and  Silberman,  (Reference  38)  obtained  the  analytical  solution 
for  a  laminar  slot  jet  heat  transfer,  and  for  a  finite  pi ate-to-nozzle  distance. 
The  final  results  are  given  graphically,  as  apparently  no  closed  form  solution 
was  obtainable.  The  Nusselt  number  for  Zn/B0s2  showed  a  small  maximum  away 
from  the  stagnation  point,  a  phenomenon  similar  to  that  found  by  Kezios 
before  for  round  impinging  jets. 
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Hrycak,  (Reference  39)  using  results  from  cold-jet  hydrodynamic  measurements, 
showed  that  a  solution  for  the  average  heat  transfer  in  the  wall-jet  region 
could  be  obtained  through  the  so-called  Colburn  analogy,  StPr^/3  =  Cf/2;  that 
particular  solution  included  also  the  stagnation  region: 

Mu  «  1.95  Pr1/3Re00-7  (D/Dp)1,23  (18) 

and  is  independent  of  the  Zn/D  term  for  Zn/0<  7.  The  above  expression  applies 
to  a  single  round  jet  of  nozzle  diameter  0  impinging  normally  at  the  center  of 
a  round  plate  of  diameter  Dp.  Other  authors  used  this  approach  before,  but 
Equation  (18)  is  based  on  a  Cf  specifically  derived  for  waM  jets,  and  also 
uses  other  parameters,  characteristic  for  axisymmetric  wall  jets  (cf.  also 
Reference  5). 

From  looking  at  the  available  analytical  solutions,  it  is  seen  that  they 
all  have  certain  common  features.  The  general  form  of  the  equation  for  the 
Nusselt  number  is 

Nu  =  C  Prm  Ren  (19) 

with  m  ranging  from  1/3  to  0.4  and  n  ranging  from  0.5  at  the  stagnation  point 
to  0.7  and  more  in  the  wall-jet  region. 

The  value  of  C  varies  from  0.763  (Reference  28)  to  1.09  (Reference  36) 
for  circular  jets,  and  Is  1.14  for  slot  jets  (Reference  31),  according  to 
Eckert's  method.  For  the  wall-jet  region,  a  reasonably  good  agreement  between 
theory  and  experiment  has  been  found  .  For  the  cases  of  even  greater  complexity, 
dimensional  analysis  appears  to  be  very  useful. 

3.2  Dimensional  Analysis 

In  many  physical  situations  involving  Impinging  jets,  the  geometry  Is  so 
complicated  that  the  governing  differential  equations  cannot  be  solved  at  all 
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analytically.  Dimensional  analysis  is  very  helpful  in  such  cases;  additionally, 
It  Is  very  handy  In  correlating  experimental  data,  by  reducing  the  total  number 
of  the  variables.  Dimensional  analysis  teaches  us  (Reference  12)  that  heat 
transfer  from  Impinging  jets  can  be  correlated  in  the  form 

Nu  •  f(Re,  Pr,  M,  r/D,  Zn/D,  Dp/D,  Dc/D,  Cn/D,  L-j/D,  L2/D...Ln/D)  (20) 

which  results  can  also  be  represented  for  restricted  ranges  of  the  parameters 
Re,  Pr,  and  M  by  the  relation 

Nu  -  C  Rea  Prb  Mc  (r/D)d  (Zn/D)e  (Dp/D)f  (Dc/D)9  (21) 

(Cn/D)h  (L1/D)1...(Ln/D)u 

where  C  is  a  constant.  The  parameters  Re,  Pr  and  M  represent  the  physical  side 
of  the  problem.  For  M  <  0.5,  the  influence  of  Mach  number  is  small,  and  can 
be  safely  neglected  (here,  formally  by  setting  c  =  o).  For  boundary  layer  type 
laminar  flow,  a  *  0.5  and  b  »  0.4  (Reference  28).  For  turbulent  flow  in  the  wall- 
jet  region,  a  =  0.7  and  higher  applies,  and  b  =  0.33  is  commonly  used.  The 
other  dimensionless  ratios  are  characteristic  lengths,  normalized  with  the 
nozzle  diameter  D  for  the  sake  of  convenience.  Local  values  of  Nu  at  stagnation 
point  (r/D  <  1)  may  be  expected  to  depend  on  Zn/D  and  Dc/D,  or  on  the  normalized 
nozzle-to- target  plate  distance  and  the  calorimeter  (or  heat  flux  gauge)  dia¬ 
meter,  but  would  not  depend  on  the  target  plate  diameter,  Dp. 

Occassionally ,  other  characteristic  lengths  are  useful,  like  the  thickness 
of  the  wall  jet  boundary  layer,  or  the  jet  spreading  pattern  parameter  b.  Perusal 
of  the  pertinent  literature  teaches  that,  in  many  cases,  information  on  the 
significance  of  certain  parameters  pertinent  to  the  problem  at  hand  may  be  lack¬ 
ing,  thus  making  the  value  of  information  available  greatly  reduced  for  various 
practical  applications.  In  particular,  for  the  average  heat  transfer  results, 
the  effect  of  the  ratio  D/Dp  is  often  not  emphasized  sufficiently  well.  On 


the  other  hand,  for  strictly  local  heat  transfer  results*  one  can  expect 
some  dependence  of  the  local  Nusselt  number  on  the  relative  size  of  the 
heat  flux  gauge  or  calorimeter,  In  the  form  of  (D_/D)9.  Likewise,  for  a  row 

v* 

of  jets,  there  will  be  Nusselt  number  dependence  on  the  nozzle- to-nozzle 
spacing,  or  (Cn/D)h. 

Dimensional  analysis  Is  Indispensable  for  a  better  understanding  of  the 
literature  on  the  subject  of  impinging  jet  heat  transfer.  Reports  and  papers 
concerned  with  experimental  results  obtained  by  various  investigators  often 
appear  to  be  hard  to  understand  and  contradictory  to  the  reader.  Dimen¬ 
sional  analysis,  properly  applied,  Is  then  of  great  help.  In  a  new  field 
like  Impinging  jet  heat  transfer,  a  certain  lack  of  generality  In  the  descrip¬ 
tion  of  Its  associated  experimental  activities  is  inevitable,  and  may  be 
considered  as  a  healthy  sign  of  growth.  As  the  interest  in  Impinging  jet 
heat  transfer  grows,  one  can  also  expect  that  additional  terms  may  be  con¬ 
sidered  to  be  of  significance  In  equations  21  and  22.  For  example,  (L-j/D)’ 
term  may  be  considered  as  representative  of  the  roughness  of  the  target 
plate,  etc. 
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SECTION  IV 


EXPERIMENTAL  HEAT  TRANSFER 
4.1  Single  Circular  Jets 

Experimental  results  obtained  from  the  study  of  heated  jets  striking  cool 
flat  surfaces  are  reported  In  References  40  to  44.  Perry  (Reference  40)  con¬ 
sidered  air  heated  to  873  K  (1112°F)  discharging  through  nozzles  16.5  and 
21.6  m  (0.65  and  0.85  In.)  In  diameter  at  velocities  up  to  91  meters  per  second 
(300  ft/sec).  The  nozzle-to-plate  spacing  was  at  least  8  nozzle  diameters.  The 
effect  of  the  angle  of  jet  Impingement  was  studied.  The  heat  flow  was  determined 
by  means  of  a  16.5  mm  (0.65-in.)  diameter  calorimeter;  this  dimension,  D  ,  was 
used  as  the  characteristic  dimension  in  both  the  Nusselt  and  Reynolds  numbers. 

The  test  results  are  shown  In  Figure  4.  They  were  correlated  for  a  normal 
Impingement  by  the  expression 

Nu  =  0.181  Re0*7  Pr0,33  (22) 

Heat  transfer  was  tested  also  at  other  angles  of  impingement.  Perry  found  a 
twofold  Increase  In  the  plate  heat  transfer,  as  the  jet  axis  of  symmetry  was 
turned  from  15°  to  90°  with  respect  to  the  plate. 

The  Influence  of  nozzle-to-plate  spacing  on  the  heat  transfer  of  an  Im¬ 
pinging  hot  air  jet  on  a  flat  plate  was  Investigated  by  Thurlow  (Reference  41). 
Values  of  Zp/D  »  10  and  Re^  up  to  60  000  were  considered.  Nozzles  2.54  and  12.7 
mm  (1  and  1/2  In.)  In  diameter  were  used  to  Impinge  hot  air  on  a  copper  plate 
with  linear  dimensions  of  610  by  152  mm  (24  by  6  In.).  The  data  were  correlated 
by 

Nu  »  C  ReJ/3  exp(-0.037  Zp/D)  (23) 

where  the  coefficient  C  was  found  to  be  equal  to  1.06  for  the  25.4  mm  diameter 
nozzle,  and  equal  to  0.33  for  the  12.7  irm  diameter  nozzle,  or  C  -  D3/^. 

The  experimental  results  of  Huang  (Reference  42)  on  turbulent  hot  air  jets 
with  temperature  ranging  from  300  to  350°F,  Indicate  a  local  value  (near  the 
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'.l.rKjf  i.H  ion  point;  ol  „  ,, 

1/3  °.87 

Nu0  =  0.023  Pr  ReQ 

(24) 

and  the  average  value 

1/3  0.87 

Nu  =  0.0180  Pr  ReD 

(25) 

for  1<  Zn/D  <  10  and  for  Reynolds  numbers  going  up  to  10^.  No  variation  of 

heat  transfer  with  Zn/o  was  detected  for  Zn/D  <  6  due  most  likely  to  the  local 

2 

heat  flow  detector  sensitive  area  being  1  in.  «  with  nozzles  ranging  from  1/4 
in.  to  1/8  in.  in  diameter,  so  that  only  average  heat  flux  was  in  fact  measured 
The  experimental  results  showed  that  the  maximum  heat  transfer  was  achieved 
with  a  4.76  mm  (3/16-in.)  diameter  nozzle  and  an  r/D  ratio  less  than  3.  It 
should  be  pointed  out  that,  in  this  investigation,  no  change  was  found  in  the 
heat-transfer  coefficient  for  Zn/t)<  6;  in  contrast,  major  changes  were  found 
in  Reference  7  (to  be  discussed  below). 

Popiel  et  al . ,  in  Reference  43,  show  local  heat  transfer  from  a  normal 
impinging  hot  gas  jet  to  a  plane  isothermal  surface,  at  low  Reynolds  numbers. 
Reasonably  good  agreement  is  found  with  Sibulkin's  solution  (Figure  5),  for 
Zn/D  < 5.  Radial  heat  transfer  distributions  were  found  to  be  qualitatively 
consistent  with  other  investigations  of  the  impinging  cold  jets  at  low  Reynolds 
numbers.  Comfort  et  al.  (Reference  44)  also  investigated  hot  jets. 

Room- temperature  jets,  impinging  on  heated  flat  surfaces,  are  reported  in 
References  7,8  and  45-46.  Reference  7  presents  results  of  tests  of  cool  air 
impinging  from  a  circular  jet  onto  an  isothermal  hot  surface  by  Gardon  and 
Cobonpue.  Nozzles  with  diameters  from  2,3  to  9.0  mm  (0.089  to  0.354  in.)  and 

nozzle  exit  velocities  up  to  sonic  were  used  In  conjunction  with  a  152.4  mm 
(6-1n.)  square  aluminum  plate.  Heat-flow  measuring  gauges  0.9  mm  (0.035  in.) 
in  diameter  were  mounted  flush  in  the  aluminum  plate  to  measure  heat-transfer 
rates.  In  a  later  publication,  (Reference  9)  Gardon  stated  that  calibration 
of  his  heat  flow  gauges  could  have  been  as  much  as  40  percent  too  high.  Provisions 
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were  made  to  Investigate  various  nozzle-to-plate  spacings  (0.25  to  40  diameters). 

Stagnation  point  heat  transfer  presented  In  Reference  7  had  the  Nusselt 
number  based  on  the  nozzle  diameter  and  the  Reynolds  number  on  the  nozzle 
diameter  and  nozzle  exit  velocity.  The  results  are  shown  in  Figures  5  and  6 
for  various  nozzle  diameters  and  nozzle-to-plate  spacings.  The  figures  show  for 
peak  heat-transfer  rates  a  tendency  to  occur  at  progressively  lower  nozzle- 
to-plate  spacings,  the  peaks  Increasing  with  Increasing  nozzle  diameter. 

Figure  6  shows  also  clearly  the  residual  effects  of  nozzle  diameter  on  heat 
transfer  for  all  Reynolds  numbers;  this  effect  Is  most  pronounced  for  Zp/D  <  10, 
however.  For  values  of  Zn/D  >  20  and  nozzle  exit  Reynolds  number  bigger  than 
14  000,  the  following  correlation  was  obtained: 

Nuq  -  13  D/Zn  (26) 

The  magnitude  of  heat  transfer  Intensity  reported  by  Gardon  and  Cobonpue  was 
substantially  less  than  that  found  by  Comfort  et  al .  (Reference  44)  for  stag¬ 
nation  point  heat  transfer,  relatively  small  Zp/D,  and  M=0.762. 

The  cooling  of  flat  circular  plates  by  air  jets.  Impinging  normally  on  the 
plates,  was  Investigated  In  Reference  30  by  Walz.  The  transient  technique  was 
used  to  determine  average  heat-transfer  coefficients  for  nozzle  Reynolds  numbers 
up  to  31  000,  nozzle-to-plate  spacings  from  1  to  8  nozzle  diameters,  nozzle 
diameters  from  3.2  to  7.6  mm  (0.125  to  0.3  In),  and  target  diameters  from  12.7 
to  25.4  mn  (0.5  to  1  In.).  Results  similar  to  the  original  results  of  Reference 
7  were  obtained  (before  announcement  of  the  possibility  of  Inaccurate  calibra¬ 
tion  of  the  heat  flux  gauges  used  In  Reference  7).  The  effect  of  the  ratio  of 
target  diameter  to  nozzle  diameter  Is  also  discussed. 

Experimental  heat-transfer  coefficients  were  presented  by  Brdllk  and  Savin 
In  Reference  36,  for  a  Jet  Impinging  on  a  heated  plate.  The  data  showed  that 
for  Zn/D  <  6.2,  the  heat  transfer  coefficient  was  practically  Independent  of 
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the  distance  from  the  jet  exit,  and  for  Zn/0>6.2,  the  distance  Zn  had  an 
appreciable  effect  on  the  heat  transfer.  Good  agreement  was  obtained  when  the 
experimental  and  analytical  results  were  compared. 


Room-temperature  turbulent  air  jets  impinging  on  a  segmented  flat  Invar 
steam-heated  plate  perpendicular  to  the  jet  axis  were  studied  by  Chamberlain 
(Reference  8).  Two  nozzle  sizes,  6.35  and  9.52  mm  (0.250  and  0.375  in.)  in 
diameter;  nozzle  exit  Reynolds  numbers  up  to  67  000;  and  no2zle-to-plate  spacings 
up  to  50  nozzle  diameters  were  investigated.  When  the  value  of  no2z1e-to- 
plate  spacings  Zn/D  was  7  or  less  (the  plate  was  located  within  the  potential- 
core  region),  local  stagnation-point  Nusselt  numbers  Nu0  were  correlated  by 

NuQ  =  1.16  Re0'447  Pr0,333  {21] 

When  the  plate  was  located  outside  the  potential -core  region,  the  data  were 
correlated  by 

NuQ  ■  0.384  Re°‘569  Pr0,333  for  Zn/D  >7  (28' 

where  Rea  is  the  Reynolds  number  based  on  the  arrival  velocity  (as  proposed 
in  Reference  7  when  multiple-jet  data  were  correlated).  Results  of  thermal 
conductivity  test  conducted  by  the  National  Bureau  of  Standards  on  a  sample  of 
invar  used  in  the  test  plate  showed  that  it  was  40  percent  higher  than  that 
available  from  the  literature  previously.  The  coefficients  in  the  two  preceding 
equations  have  been  modified  accordingly. 

For  locations  along  the  plate  away  from  the  stagnation  point,  the  correlation 


~=  exp.  (-  1.56  (r/Znr  /D)  (29) 

no 

was  suggested;  here  hr  is  the  local  coefficient  at  a  distance  r  from  the  stagnation 


point. 


Heat-transfer  coefficients  for  a  jet  impinging  on  a  plate  are  reported  In 
Reference  45,  for  jet  exit  velocities  of  61  to  213  meters  per  second  (200  to 
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700  ft/sec).  Calorimeters  were  used  tc  obtain  the  heat-transfer  distributions 
on  the  plate.  Turbulence  characteristics  of  the  applicable  free  jets  were 
used  as  a  basis  for  correction  factors  to  be  applied  to  computed  laminar  stagnation 
point  heating. 

In  particular,  near  the  stagnation  point,  the  heat  transfer  may  be  computed 
from  the  laminar  heat  transfer  that  would  take  place  on  a  surface  having  the 
same  pressure  distribution  as  that  on  the  impingement  plate  and  using  a  correction 
factor,  while  farther  away  from  the  stagnation  point,  the  heat  transfer  may  be 
computed  using  the  technique  of  Reference  46,  and  very  far  from  the  stagnation 
point,  the  relation 

Nu  =  0.12  Rey®  (30) 

should  be  used,  where  the  Reynolds  number  is  based  on  the  distance  at  which 

the  velocity  is  one-half  the  maximum  free-jet  velocity  in  the  plane  of  impingement. 

Interesting  measurements  of  heat  transfer,  from  naturally  occuring  and 
turbulized  jets,  have  been  reported  by  Dyban  and  Mazur  in  Reference  47.  They 
show  clearly  the  enhancing  effects  of  turbulence  on  stagnation  point  heat 
transfer  in  Figure  7.  Heat  transfer  boosts  of  the  order  of  100  percent  could 
be  produced  by  turbulizing  the  impinging  jet,  that  is  by  increasing  artificially 
the  fluctuating  velocity  components  in  the  jet  stream.  Additionally,  Dyban  et 
al.in  Reference  48  show  the  effects  of  turbulizing  on  the  term  a*,  which  is 
needed  in  analytical  heat  transfer  calculations  of  impinging  jets  (Figure  8). 

In  Figure  8,  comparisons  with  results  from  Reference  5,  11  and  48  are  made. 


It  is  seen  that  a*  may  be  approximated  by 

a*  =  (Zn/D)0J6,l  <Zn/D  <6.2 

(31) 

and 

a*  »  82.6  (Zn/D)_1'75,Zn/D  >6.2 

(32) 

Comparing  the  results  of  Poplel,  Gardon  and  Cobcnpue,  and  those  of  Dyban  et  al . 
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In  Figures  5-8,  it  is  seen  that  turbulence  effects  enhance  heat  transfer  by 
nearly  100  percent  at  the  stagnation  point,  in  case  of  impinging  round  jets, 
in  proportion  to  (Z  /D)®*^  for  Zn/D  less  than  6  plus.  The  dimensionless 
velocity  gradient  at  the  stagnation  point  is  also,  at  first,  increasing  with 
Zn/D,  until  Zn/D  =  6.2  Is  reached;  then  a  gradual  decrease  in  a*  is  observed. 

The  effects  of  turbulence  Influence  apparently  the  rate  of  fluid  mixing  in 
the  "mixing  cone"  of  the  jet,  and  also  the  length  of  the  potential  core.  These 
facts  influence  pressure  gradients  at  the  target  plate  in  the  direction  parallel 
to  the  plate,  which  explains  the  changes  in  a*  shown  in  Figure  8.  In  Figure 
6,  it  is  seen  that  many  nozzles  achieve  a  high  degree  of  turbulence  in  the 
jets  they  discharge,  without  any  apparent  use  of  turbulizers.  This  may  be  con¬ 
cluded  by  comparing  the  Figure  6  heat-transfer  data  with  those  in  Figure  7, 
where  turbulizers  were  purposely  used. 

4.2  Arrays  of  Circular  Jets 

Reference  1  presents  results  of  an  experimental  investigation  of  heat 
transfer  by  Freidman  and  Mueller,  for  an  array  of  circular  air  jets  Impinging 
on  a  heated  plate.  Effects  of  hole  size  and  spacing,  distance  between  the 
array  of  nozzles  and  the  heated  plate,  and  air  velocity  on  the  heat-transfer 
rates  were  studied.  Hole  sizes  from  6.34  to  19.05  mm  (0.250  to  0.750  In.)  and 
spaclngs  between  the  nozzles  and  the  heated  plate  from  57.2  to  158.8  mm  (2.25 
to  6.25  In.)  were  considered.  Results  indicated  that  spaclngs  between  holes 
of  from  4  to  6  hole  diameters  (giving  about  5  to  2  percent  free-flow  area) 
yielded  the  best  heat  transfer  results. 

In  addition  to  the  study  of  single-jet  impingement  discussed  above,  Gardon 
and  Cobonpue  also  considered  multiple  arrays  of  jets;  the  results  are  reported 
In  Reference  7,  where  heat-transfer  effects  of  nozzles  with  diameters  from 
1.59  to  12.7  mm  (1/16  to  1/2  in.),  arranged  in  square  arrays,  were  studied. 
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A  5x5  nozzle  array  on  50.8  nm  (2-1n.)  centers  and  a  7x7  nozzle  array  on  30.5 
mm  (1.2-in.)  centers  were  used.  Airflow  rates  from  50  to  800  1b/(hr  ft^)  were 
considered.  Nozzle  exit  velocities  were  kept  low  enough  so  that  compressibility 
effects  could  be  neglected. 

When  attempts  to  correlate  the  data  based  on  nozzle  exit  velocity  failed, 
Gardon  and  Cobonpue  (Reference  7)  decided  to  correlate  by  using  the  velocity 
of  arrival.  Space-averaged  heat-transfer  data  were  successfully  correlated 
when  the  center-to-center  hole  spacing  was  used  as  the  characteristic  dimension 
in  the  Nusselt  and  Reynolds  numbers,  and  the  velocity  of  arrival  was  used  In 
the  Reynolds  number.  The  velocity  of  arrival  ua  for  nozzle-to-plate  spaclngs 
greater  than  8  was  evaluated  as 

ua  =  6.63  uQc  D/Zn  for  Zn/D  >  8.  (33) 

The  average  Nusselt  number  was  found  to  be 

Nu  *  0.286  Re°‘625  (34) 

a 

The  equation  for  Nu  above  holds  down  to  Zn/D  >  1. 

The  effect  of  nozzle-to-plate  separation  distance  on  the  local  Nusselt 
number  was  also  investigated.  It  was  found  that  the  shapes  of  the  curves, 
notably  different,  depend  on  Rea  as  well  as  on  Zn/Cn  and,  when  Zn/D  <  8,  on 
Zn/D  also.  In  Reference  42,  the  maximum  heat  transfer  for  a  single  jet  was 
found  for  a  4.76  mm  ( 3/1 6-1 n . )  diameter  nozzle  and  an  r/D  <  3;  Huang's  investi¬ 
gation  of  multiple  jets  was  limited  to  the  4.76  mm  (3/16-1n.)  diameter  nozzle, 
spaced  at  Cp  =  15.88  ran  (5/8  in.)  center-to-center  in  a  single  row  and  the 
row-to-row  distances  of  38.1,  76.2  and  152.4  mm  (1.5,  3,  and  6  in.). 

Huang  proposed  correlations  that  came  under  a  heavy  criticism  during  the 
discussion  of  his  paper,  since  average  heat  transfer  coefficients  from  his 


arrays  of  round  jets  seemed  to  be  higher  than  those  of  the  equivalent  number 
of  single  jets.  During  discussion  Huang  stated  that  his  heat  transfer  results 
from  arrays  of  jets  wore  carefully  determined.  Comparison  of  his  data  is 
shown  in  Figure  10. 

A  triangular  array  of  circular  jets,  impinging  on  a  flat  surface,  was 
investigated  by  Ott  in  Reference  50.  Only  a  limited  range  of  Reynolds  numbers 
was  considered  and  only  Zn/D  and  Re  were  variable.  An  average  value  of  the 
heat-transfer  coefficient  was  therefore  correlated  with  only  Re  and  Zn/D  as 
parameters.  No  variation  in  free-flow  area  was  included  in  this  study.  Ott's 
results  are  shown  in  Figure  10. 

Hilgeroth  (Reference  51),  discovered  that  jet  interaction  and  spent  air 
effects  had  a  substantial  effect  on  heat  transfer.  He  found  also  that, 
for  Cn/D,  Z^/D,  and  u  held  constant,  the  heat-transfer  coefficients  increased 
as  the  hole  diameter  increased.  A  25  percent  decrease  in  heat-transfer  coefficient 
between  Zp/D  of  2  and  6  was  observed.  Maximum  values  of  h  were  obtained  for 
Zn/D  of  6.3,  for  a  free-flow  area  of  1.5  percent.  The  Reynolds  number  exponent 
was  found  to  be  a  function  of  hole-spacing-to-hole-diameter  ratio,  increasing, 
as  Cn/D  increased  from  3.5  to  12,5.  A  square  array  of  jets  was  found  to  be 
superior  to  an  equilateral  triangular  array  from  the  same  jet  velocity  and 
open-area  ratio.  Comparison  of  Hilgeroth's  results  with  those  of  others  is 
also  shown  in  Figure  10. 

Reference  52  reported  the  investigation  of  a  series  of  configurations  by 
Kercher.  Figure  10  compares  Kercher’s  data  and  the  results  of  Gardon  and 
Cobonpue  (Reference  7).  The  differences  between  the  results  may  be  attributed 
to  the  influence  on  boundary- 1 ayer  phenomenon  of  the  effects  of  jet  impingement 
from  perforated  plates  and  from  long  throat  nozzles.  Data  indicated  that  the 
average  surface  heat-transfer  coefficients  increased  as  Z(1/D  was  increased 
from  1  to  about  5. 
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The  test  results  are  sumnarized  in  Reference  52  as  follows:  Among  the 
highlights  we  find  that  heat-transfer  coefficients  increase  with  increasing 
open  area,  heat  transfer  is  dominated  by  the  hole-diameter  Reynolds  number  Re^ 
and  the  hole-spaclng-to-hole-diameter  ratio;  within  the  range  tested,  increasing 
Zn/D  increases  heat  transfer  without  crossflow,  but  decreases  heat  transfer 
with  cross  flow;  increasing  cross  flow  decreases  heat  transfer  performance; 
and  decreasing  hole  diameter  with  Increasing  number  of  holes,  everything  else 
being  equal,  improves  heat-transfer  performance. 

An  experimental  study  of  the  effects  of  cross  flow  on  the  heat-transfer 
characteristics  of  single  rows  of  air  jets  impinging  on  plane  surfaces  is 
reported  by  Metzger  and  Korstad  in  Reference  53.  Rows  of  10  circular  jets, 
each  2.54  mn  (0.1  in.)  in  diameter,  with  center-to-center  spacing  ratios  Cn/D 
from  2.5  to  5,  were  considered.  Nozzle-to-plate  spacings  Zn/D  from  2  to  6.7 
and  ratios  of  cross  flow  jet  flow  M*  from  1  to  3  were  studied.  The  average 
Stanton  number  was  defined  as 


St  *  St  (Ree$,  Cp/D,  2n/D,  1/D)  (35) 

where  St  Is  h/(GCp)  evaluated  over  the  cooled  surface  of  half-length  1  and 

Re  is  based  on  the  width  of  a  slot  whose  area  is  equal  to  that  of  the  10 
es 

holes  used,  and  where  G  =  Vp,  the  main  flow  rate  per  unit  area,  and  where  Bg$  = 
equivalent  slot  width,  or 


Rees  =  26  V* 


(36) 


An  extra  parameter,  M*  =  ,nc/mj »  should  correlate  data  with  cross  flow.  As  M* 
Increases,  the  general  trend  shows  a  decrease  in  St/St^  _  q  In  the  region 
upstream  from  the  nozzle  exit  plane  and  an  increase  In  the  downstream  region. 
For  the  ranges  of  variables  covered  in  this  investigation,  the  cross  flow  heat- 


transfer  results  were  correlated  by 
St  =  0.822  M*-0-049  ReD"0,338 


(37) 
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The  significant  length  on  which  Nu  and  Re  are  based  may  be  the  nozzle 
diameter  or  the  radius  or  the  target  plate  diameter  or  radius.  For  arrays  of 
jets,  the  nozzle-to-nozzle  spacing  Cn  is  often  chosen.  Most  impinging  jet  re¬ 
sults  are  stated  in  terms  of  ReD<  Equation  (34)  (for  Zn/D<8)  may  be  con¬ 
verted  into  expressions  involving  ReD  by  multiplying  each  side  by  D/Cn ,  or 

(D/Cn)  hCn/k  =  0.286  (UocCnp/y)°‘625(D/Cn) 
with  the  result 

hD/k  =  Nu  “  0.286  (Cn/D)'0,375  ReD0,625  (38) 

Equation  (38)  allows  us  to  introduce  here  explicitly  the  effects  of 
nozzle-to-nozzle  spacing  and  of  the  diameter  size  in  the  form  of  the  ratio 
(Cn/D)  to  some  negative  power.  It  indicates  that,  as  Cn  increases,  there  are 
fewer  nozzles  per  unit  area,  and  average  heat  transfer  is  reduced.  Likewise, 
with  an  increasing  nozzle  diameter,  0,  average  heat  transfer  improves  as  the 
nozzle  area  increases  in  proportion  to  the  total  area  exposed  to  impinging  jets. 

The  influence  of  the  size  of  the  nozzle  array,  D^,  can  be  also  introduced  in 
Equation  (38)  by  multiplying  each  side  with  (D/Deq)n,  n  being  obtained  experi¬ 
mentally,  where  Dgq  =  2  (area  of  nozzle  array/*)*5. 


In  this  way,  for  example,  Kercher’s  result 


Nu-  =  0.324  Re°‘658 

(39) 

°n  Ln 

can  be  transformed  into 

Nu0  =  0.324  (Cn/D)"°*342  Re0-658 

(40) 

to  be  compared  with  the  results  of  other  investigators  based  on  the  nozzle 
diameter  D,  in  Figures  9  and  10.  Equations  (38)  and  (40)  allow  us  to  Intro¬ 
duce  the  effect  of  Cn  and  the  nozzle  diameter  D  explicitely  Into  the  correlating 
equation.  The  effect  of  (D/DeC( )n ,  not  shown  specifically  in  above  examples,  is 
useful  in  comparison  of  nozzle  arrays  of  various  sizes.  Its  function  is  similar 
to  that  of  the  term  (D/Dp)"'’23  in  Equation  (18). 
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It  is  sometimes  very  convenient  to  use  the  analogy  between  heat  and  mass 
transfer  to  overcome  some  of  the  difficulties,  inherent  In  getting  highly 
accurate  measurements  of  local  heat  fluxes  due  to  multiple  jet  impingement. 

Thus,  Koopman  and  Sparrow  (Reference  54)  have  applied  the  naphtalene  sublima¬ 
tion  technique  in  conjunction  with  the  semi -automated  data  acquisition,  to 
determine  the  local  and  average  mass  transfer  coefficients  resulting  from  im¬ 
pingement  of  a  row  of  jets  onto  a  flat  surface.  The  local  Sherwood  numbers  thus 
obtained  could  then  be  used  to  Identify  the  locations  where  the  impinging  jets 
would  provide  the  most  effective  heating  or  cooling.  It  has  been  observed,  in 
particular,  that  because  of  the  collisions  of  the  wall  jets  originating  from  the 
adjacent  impinging  jets,  "the  transfer  coefficients  at  positions  mid-way  between 
the  jets  may  take  on  relatively  large  values".  The  effect  of  such  collisions  is 
enhanced  at  high  Reynolds  numbers,  small  jet  orifice  spaclngs,  and  small  oriflce- 
to-target  distances.  A  two-dimensional  wall  jet  is  produced  at  sufficiently 
large  downstream  distances  from  the  row  of  jets. 

In  contrast  to  single  jets,  where  the  results  depend  only  on  the  radial  dis¬ 
tance  to  the  stagnation  point,  for  multiple  jets  transfer  coefficients  depend 
differently  on  the  spannwlse  and  streamwise  relative  distances.  Under  certain 
circumstances,  the  multiple  jet  transfer  coefficients  were  found  to  exceed  those 
for  the  single  jet. 

On  the  other  hand,  when  the  single  jet  results  from  Reference  54  were  com¬ 
pared  with  those  of  Gardon  and  Akflrat  (Reference  9),  the  agreement  was  found 
to  be  satisfactory  for  conditions  corresponding  to  common  values  of  the  character¬ 
istics  parameters. 

In  spite  of  proven  high  precision  of  mass  transfer  results  at  lower  Reynolds 
numbers,  It  remains  doubtful  If  such  methods  may  be  considered  entirely  satis¬ 
factory  for  really  high  Reynolds  numbers  and  for  determination  of  the  maximum 
possible  heat  fluxes  in  cases  of  complex  geometry  existing  in  some  applications  of 
impinging  jets.  Such  applications  are  described  in  some  detail  below. 
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4.3  Single  Slot  Jets 

An  experimental  study  of  the  heat-transfer  characteristics  of  slot  jets 
impinging  normally  on  a  flat  target  was  performed  by  Metzger  (Reference  31). 
Nozzles  with  lengths  of  19.0  mm  (0.75  in.)  and  widths  of  0.25,  0.51,  1.02, 
1.53,  and  2.04  mm  (0.01,  0.02,  0.04,  0.06,  and  0.08  in.)  were  studied  as  they 
Impinged  on  two  copper  targets  of  half-lengths  12.7  and  6.35  mm  (0.50  and  0.25 
In.)  mounted  In  blocks  of  balsa  wood.  Reynolds  numbers  up  to  10  000  (based  on 
nozzle  exit  velocity  and  nozzle  hydraulic  diameter)  were  considered.  Initial 
tests,  with  all  variables  except  nozzle-to-target  spacing  held  constant,  show¬ 
ed  that  maximum  average  heat  transfer  occurred  at  a  value  of  Zn/BQ  =  8.  Since 
no  significant  variation  was  noted  as  Zn/BQ  was  varied  from  7  to  10,  It  was 

decided  that  all  other  tests  would  be  made  for  ZVEL  *  8.  This  value  of  8  Is 

n  o 

comparable  with  the  value  of  Zn/D  of  between  6  and  7  for  a  single  circular 
jet. 


The  correlation  equation  for  the  average  heat  transfer  was  found  to  be 

St  (ReDh)°‘434  Pr°'63  *  0.74  (0p/Bo)‘0,434  (41 

for  7  <  Zn/B0  <  10  and  3  <  Dp/BQ  <  50. 

Experimental  wall  jet  heat-transfer  data  that  were  obtained  by  Myers  et 
al.  using  a  12.7  crm  (0.5-1n.)  wide  slot  and  slot  Reynolds  numbers  from  16  600 
to  38  100  are  reported  ;n  Reference  32.  Data  were  taken  for  values  of  r/BQ 
from  20  to  186  for  five  different  unheated  starting  lengths  l  (r  Is  the  dis¬ 
tance  along  the  wall).  The  final  results  are  represented  by  the  formula 
good  for  correlation  of  the  data  for  r/B0>45,  namely. 


St  Re 


0.2 

B„ 
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(42) 


Heat-transfer  characteristics  of  air  jets  issuing  from  slots  1.59,  3.17, 
and  6.35  mm  (1/16,  1/8,  and  1/4  in.)  wide  and  152.4  mm  (6  in.)  long  and  impinging 
normally  on  an  aluminum  152.4  mm  (6  in.)  square  plate  are  reported  by  Gat  don 
and  Akfirat,  Reference  9.  Both  average  and  local  heat-transfer  rates  were 
obtained,  with  the  help  of  a  0. 9-millimeter  (0.035  in.)  diameter  heat  transfer 
gauge.  Maximum  stagnation  heat  transfer  was  found  to  be  a  function  of  the 
nozzle  exit  Reynolds  number  and  occurred  between  values  of  Zn/Bo  from  7  to  10. 

For  turbulent  jets  ( Re>  2000 )  and  nozzle-to-plate  spacings  greater  than  14 
slot  widths,  the  stagnation-poing  heat-transfer  coefficients  were  correlated 
within  5  percent  by  _q  ^ 

Nu„  =  1-2  Re“  SS  <Z„/80)'  '  (43) 

°0 

ror  ReD  up  to  50  000  and  Zn/B„  up  to  60.  Nozzle  size  was  found  to  Influence 
Bo  0 

the  turbulence  levels  of  the  jets;  for  ReBo  =  11  000,  the  Initial  turbulence 
levels  were  found  to  be  0.6,  2.5,  and  7.5  percent  for  the  1.59-,  3.17-,  and 
6.35-mm  (1/16-,  1/8-,  and  1/4-in.)  wide  slots,  respectively.  The  observed 
variation  of  the  local  heat-transfer  coefficients  along  the  plate  has  been 
dependent  on  Zn/BQ.  Thus,  for  nozzle-to-plate  spacings  Zn/B0  greater  than  14 
slot  widths,  the  variation  of  heat-transfer  coefficients  has  a  bell  shape, 
while  for  14<Zn/B0<8,  the  bell  shape  is  modified  slightly  by  an  abrupt  change 
in  slope  in  the  vicinity  of  r/B0  *  4,  On  the  other  hand,  as  Zn/Bo  Is  reduced 
below  8,  two  "humps"  begin  to  form  at  about  r/B0  =  7;  and  for  Zn/B0<  6,  they 

become  well-defined  secondary  peaks  in  the  heat-transfer  rate.  These  secondary 
peaks  are  tentatively  explained  as  due  to  a  transition  from  laminar  to  turbulent 
boundary  layers,  and  as  Zn/B0  is  reduced  below  1/2,  the  impinging  jet  becomes 
essentially  a  "wall  jet".  Heat-transfer  coefficients  rise  sharply  with  Increasing 
velocities  in  the  gap  between  the  nozzle  exit  and  the  plate;  the  local  heat- 
tranfer  coefficients  show  two  peaks  away  from  the  stagnation  line. 
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Avnr.iyn  ho.it- transfer  coefficients  for  the  three  slot  widths  considered 
are  correlated  against  the  Reynolds  number  based  on  the  arrival  velocity  by 
an  expression  of  the  form. 

Nu  *  0.36  Rea’^ 
for  Rea>10^,  and  Zn/B0>8, 

The  heat-transfer  characteristics  for  a  slot  jet  impinging  on  a  flat  copper 
surface  for  slot  widths  of  6.35,  12.7  and  19.0  187.5  mm  (1/4,  1/2,  and  3/4 
in.),  nozzle-to-pl ate  spacings  from  2  to  16  slot  Reynolds  numbers  from  4600 
to  102  000  were  studied  by  Cadek  (Reference  33).  Local  heat-transfer  rates 
were  measured  by  use  of  a  miniature  circular  foil  heat-flow  sensor  with  a  0.91 
mm  (0.036-in.)  diameter  sensing  area.  Values  of  r/B0  from  0  to  36  were  considered. 
Maximum  heat  transfer  was  found  for  Zn/B0  s  8,  in  agreement  with  Gardon  and 
Akfirat's  data.  Cadek  concluded  from  his  study  that  for  2<Zn/Bo<il*  9°°d 
agreement  between  theory  and  experiment  was  obtained  for  both  the  stagnation 
region  and  the  wall  jet  region. 

Reference  55  presents  experimental  heat-transfer  data  obtained  by  Cartwright 
and  Russel,  from  an  impinging  jet  for  nozzle-to-plate  spacings  from  8  to  47 
slot  widths,  nozzle  exit  Reynolds  numbers  from  25  000  to  110  000,  and  distances 
along  the  wall  jet  from  0  to  132  slot  widths.  At  low  slot  Reynolds  numbers, 
the  maximum  value  of  the  heat-transfer  coefficient  occurred  at  the  stagnation 
point,  and  then  decreased  monotonical ly  with  increasing  distance  from  the 
stagnation  point,  while  at  larger  Reynolds  numbers,  the  maximum  heat-transfer 
coefficient  occurred  some  distance  from  the  stagnation  point,  in  a  way  similar 
to  axisymmetric  impinging  jet  heat-tranfer  results. 

Analytical  attempts  to  predict  the  heat  transfer  appeared  to  have  been 
successful  only  In  the  wall  jet  region  beyond  r/B0  *  36;  failure  of  predictions 
in  the  impingement  region  was  explained  by  the  high  level  of  turbulence  existing 
in  the  free  jet  before  impingement. 
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An  analytical  method  for  predicting  the  local  Nusselt  number  In  the 
stagnation  region  of  a  slot  jet  was  developed  by  Andreyev  et  al.  in  Reference 
56.  The  prediction  was  eventually  modified,  by  use  of  some  experimental  data 
so  that  an  empirical  correlation  resulted,  where  e  Is  the  degree  of  turbulence 


Intensity  for  1  <  Zn/Bc  <  6.5: 


Nu  =  0.48  - Vr  [  1 

<W 


0.116 


(r/Bj‘ 


W 


rl  (1  +  0.01 5e) 


The  expressions  for  6.5  <  Zn/B0  <  12  and  Z^/B0  >  12  are  of  the  same  form,  with 
some  changes  in  the  values  of  the  numerical  constants  occurring  In  Equation  45. 


4.4  Arrays  of  Slot  Jets 

The  tests  reported  In  Reference  1  by  Freldman  and  Mueller  included  infor¬ 
mation  on  plates  with  a  multitude  of  slots.  It  was  found  that  the  slotted  plates 
were  generally  comparable  in  performance  to  the  round  nozzle  equipped  plates, 
with  some  exceptions. 

Gardon  and  Akflrat  In  Reference  9  also  report  the  results  of  an  investiga¬ 
tion  of  Impingement  from  an  array  of  slot  jets.  Data  were  obtained  for  two 
combinations  of  slot  Jets,  namely,  three  slots  at  50.8  mm  (2-1 n . )  centers,  and 
two  slots  at  101.6  mm  (4-1n.)  centers.  It  was  found  that,  at  the  smallest 
slot-to-plate  spacing,  the  Identity  of  each  jet  was  preserved  and  the  peak 
heat-transfer  rates  of  the  jets  differed  only  slightly  from  those  corresponding 
to  single  jets.  Interaction  between  the  jets  produced  secondary  peaks  midway 
between  the  points  of  Impingement.  As  the  length  of  the  jet  was  Increased, 
Interaction  occurred  before  Impingement.  The  peaks  did  not  necessarily  lie 
below  the  centers  of  the  jets.  At  the  largest  spacing,  the  jets  behaved  al¬ 
most  like  a  large  single  jet. 
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The  data  were  correlated  on  the  basis  of  a  Reynolds  number  based  on  the 
velocity  of  arrival  defined  here  In  Equation  3  with  the  average  Nusselt  manber 
defined  as 

Nu  •  F  Cn/k  (46) 

In  the  form 

Nu  *  0.36  Re?>62  (47) 

d 

which  expression  should  not  be  extrapolated  to  Cn/BQ  <  16. 

Schuh  and  Pettersson  report  heat-transfer  data  In  Reference  57  for  Impinge¬ 
ment  of  air  through  arrays  of  slots.  Slot  widths  of  1  and  5  nm  (0.0394  and  0.197 
in.),  Jet  spaclngs  of  25,  40,  60  and  100  mm  (0.984,  1.57,  2.36,  and  3.94  In.), 
and  slot-to-plate  spaclngs  of  2,  4,  6  and  8  slot  widths  were  Investigated.  Slot 
exit  velocities  from  4  to  100  meters  per  second  (13.12  to  328  ft/sec)  were  con¬ 
sidered.  For  12  000  <  Re  ±  100  000,  5  <  Cn/BQ  <_  100,  and  Zn/B0  •  4,  It  was  found 

that  the  data  could  be  correlated  by 
-  -0.327 

St  -  0.461  (A  Reg0, 402  (48) 

o  o 

within  10  percent.  Also,  In  Reference  57,  It  was  found  for  rows  of  jets  that 
lower  heat-transfer  was  achieved  for  the  round  jets  than  for  the  arrays  of  slot 
Jets  with  wide  row  spaclngs.  Little  difference  was  found  for  arrays  with 
narrow  row  spaclngs. 

It  was  also  found  that  a  superimposed  wall-parallel  flow,  of  up  to  60 
percent  of  the  jet  flow,  could  be  Imposed  with  no  reduction  In  heat  transfer. 

For  similar  conditions,  but  with  rows  of  round  Jets  Instead  of  arrays  of  slot 
Jets,  some  reduction  In  heat-transfer  accompanied  the  60  percent  superimposed 
wall  parallel  flow. 
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In  conclusion  of  discussion  related  to  the  slot  jet  experiments,  sub¬ 
stantial  experimental  difficulties  In  obtaining  reproducible  slot  jets  should 
be  pointed  out.  While  actual  slot  jets  are  all,  by  necessity,  of  finite  length, 
great  care  must  be  taken  to  make  sure  that  the  end  effects  are  eliminated. 

For  that  reason,  the  comparison  of  Individual  slot  jet  data  must  be  made  with 
some  caution. 

Heat  transfer  from  Impinging  jets  literature  review,  carried  out  In  1968- 
72,  has  been  prepared  by  Llvlngood  and  Hrycak  (Reference  58).  This  informa¬ 
tion,  In  an  updated  version,  has  been  used  In  Sections  III  and  IV  of  this  re¬ 
port.  Another,  more  recent  literature  related  review  has  been  made  by  Martin 
(Reference  59).  It  focusses  its  attention  to  a  large  extent  on  the  results 
of  German  workers  In  the  present  area  of  Interest,  and  deals  with  heat  transfer 
optimization  problems  of  arrays  of  jets  Impinging  on  flat  plates.  Since  heat 
transfer  from  arrays  of  jets  has  not  yet  been  described  In  a  definitive  way, 
additional  work  Is  Indicated  before  heat  transfer  from  impinging  jets  can  be 
optimized  In  a  physically  meaningful  fashion.  Still,  such  efforts  are  very 
timely.  Figure  10  has  been  adapted  from  Reference  60.  In  Reference  60,  dimen¬ 
sions  of  characteristic  parameters  have  been  kept  down  to  what  may  be  expected 
In  Impingement  cooling  of  a  typical  gas  turbine  blade.  The  results  may  be 
quantitatively  different  from  other  cooling  applications,  with  substantially 
larger  diameter  nozzles,  as  Is  positively  shown  In  Figure  10. 

4.5  Cooling  of  Areas  with  a  Curvature 

Cooling  of  areas  with  a  curvature  is  an  Important  application  of  rows  and 
arrays  of  Impinging  Jets,  primarily  to  cool  gas  turbine  blades  Internally. 

Such  areas  may  be  approximated  by  concave  semi  cylindrical  surfaces.  Among  the 
Investigators  here,  we  may  distinguish  Metzger  et  al . ,  Chupp  et  al.,  Jenkins 
and  Metzger,  Oyban  and  Mazur,  Tabakoff  and  Clevenger,  and  Burggraf 
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(References  61-66  respectively).  Very  helpful  Is  here  also  the  monograph 
by  Shvets  and  Dyban,  with  an  extensive  bibliography  (Reference  67). 

In  cooling  of  areas  with  a  curvature,  several  unique  effects  must  be 
considered.  The  most  significant  effect  Is,  perhaps,  that  of  the  severely 
constricted  space  between  the  nozzle  exits  and  the  Impingement  plate  with  a 
curvature.  It  becomes  very  significant  when  the  nozzles  are  formed  by  a  row 
of  round  holes,  In  a  tightly  fitting  plenum  chamber.  Because  of  a  restricted 
space  In  the  case  of  low  Zn/D  values,  a  true  wall-jet  cannot  develop  In  all 
cases,  and  the  stagnation  area  (where  Re0*®  term  Is  still  applicable)  may  be 
somewhat  extended.  The  restricted  free  space,  and  the  fact  that  the  ends  of 
the  heat  transfer  surface  must  be  blocked  off,  to  make  flow  still  two-dimen¬ 
sional,  Imply  for  low  Zn/d  values  the  possibility  of  formation  of  a  channel- 
like  flow,  as  suggested  In  Reference  64. 

For  heat  transfer  at  the  stagnation  point,  Hrycak  (Reference  39)  developed 
an  expression 

Nuq  -  1.6  Re[J*5  Pr0,39  (Dc/D)'0*42  (Zn/D)'0,22  (Cn/D)"0*16  (49) 

that  shows  expllcltely  the  effect  of  the  nozzle  diameter,  D,  of  Zn/D  and  nozzle- 
to-nozzle  spacing,  Cn< 

The  expression  Is  similar  to  a  corresponding  flat  plate  formula  that  Is 
also  discussed  in  Reference  39.  For  the  average  heat  transfer,  Reference  39 
reports  the  expression 

Nu  «  0.72  Pr0*33  ReD0,63  (Cn/0)“0,16  (D/Dp)0,402  (50) 

Equation  (46)  Is  compared  In  Figures  9  and  10  with  arrays-of-nozzles  results  by 
Kercher,  and  Gardon  and  Cobonpue.  It  seems  to  fall  almost  exactly  Into  the 
middle  between  the  results  of  these  two  Investigators.  Equation  (50)  describes 
heat  transfer  on  a  semi  cylindrical  plate  10  In.  long  and  5  In.  in  diameter; 
comparison  with  the  results  of  other  Investigators  working  with  a  similar 
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geometry  is  shown  In  Figure  11,  while  local  distribution  of  heat  transfer 
coefficients  Is  shown  in  Figure  12.  In  both  cases,  10  by  5  in.  semi  cylindrical 
plate  results  compare  favorably  with  the  results  of  others.  Also,  heat  transfer 
on  that  10  by  5  in.  semi  cylindrical  plate  mentioned  above  has  been  found  similar 
to  heat  transfer  on  a  plate  scaled  down  ten  times.  Performance  of  that  small- 
scale  configuration  plate  has  been  compared  with  the  results  of  several  curvilinear 
geometries  In  Figure  13.  Again  the  comparison  is  quite  favorable.  This  figure 
has  been  adapted  from  an  Illustration  prepared  by  Llvingood  and  Gauntner  In 
Reference  68,  while  Figure  12  was  adopted  from  Reference  69.  In  Figures  11-13 
there  occur  references  that  have  not  been  listed  directly  In  the  present  report. 
These  additional  references  may  be  found  in  References  68  and  69. 

Impinging  jet  heat  transfer  from  jet  rows,  and  jet  arrays,  to  both  flat  and 
curved  surfaces  depends  ultimately  on  the  performance  of  the  Individual  jets,  form¬ 
ing  the  most  essential  part  of  tne  sometimes  very  complex  heat  exchanging 
systen.  Fundamental  here  is  still  the  Information  on  a  single,  Impinging, 
axlsymmetrlc  jet.  A  lot  of  Important  Information  on  behavior  of  a  single  Impinging 
jet  (applied  to  mass  transfer),  where  the  effects  of  surface  roughness  have 
been  studied  systematically  for  the  first  time,  for  example,  has  been  collected 
by  Schrader.  The  behavior  of  a  single  jet  In  cross-flow  has  been  studied  by 
Bouchez  and  Goldstein.  A  systematic  study  of  the  impinging  jet  behavior  has 
been  carried  out  in  great  detail  by  Poplel  (References  70-72).  Additional  In¬ 
sights  Into  behavior  of  rows  of  jets  Impinging  on  semi  cylindrical  surfaces  may 
be  seen  In  a  recent  publication  by  Hrycak  (Reference  73),  where  also  some  other 
fundamental  aspects  of  Impinging  Jet  heat  transfer  are  covered. 

More  Information  on  the  Interaction  between  Individual  jets  In  an  array 
appears  to  be  very  desirable.  Details  of  how  a  centrally  located  jet  interacts 
with  Its  eight  neighbors  In  a  three-by-three  square  jet  array  would  be  probably 
the  logical  first  step  towards  this  kind  of  Investigation. 
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SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  General  Considerations,  Individual  Circular  and  Slot  Jets 

Experimental  investigation  of  heat  transfer  from  impinging  jets  has  a  much 
shorter  history  than  the  experimental  Investigation  of  free  jet  flow  patterns. 
Consequently,  agreement  among  the  results  of  various  investigators  on  a  number 
of  important  aspects  of  heat  transfer  from  impinging  jets  is  still  not  as  good 
as  could  be  desired. 

Here  we  can  distinguish  three  areas  where  experimentation  has  been  carried 
out:  on  impinging  round  jets,  impinging  slot  jets,  and  arrays  of  jets.  So  far, 
the  theory  has  been  able  to  calculate  mainly  heat  transfer  at  the  stagnation 
point  under  some  simplifying  assumptions  (cf.  References  2,30,38,  73)  and  several 
reasonably  successful  attempts  have  been  made  to  do  likewise  for  the  wall  jet 
region,  (References  14,33).  Also,  measured  heat-transfer  coefficients  in  the 
stagnation  region  exceed  those  determined  by  use  of  the  several  analytical 
procedures  discussed  herein.  The  analyses  are  affected  by  a  lack  of  knowledge 
of  the  jet  turbulence  characteristics.  Further  Investigations  leading  to  a 
more  precise  evaluation  of  the  jet  turbulence  characteristics  appear  to  be 
required  for  more  satisfactory  calculations.  Additional  Investigations  of  single, 
turbulent,  room- temperature  air  jets,  impinging  on  heated  flat  surfaces,  would 
be  useful . 

The  main  contributions  of  theory  at  the  present  time  seem  to  be  the  following: 
for  efficient  heat  transfer,  for  both  circular  and  slot  jets,  the  Impinging 
surface  has  to  be  at  the  end  of  the  potential  cone,  where  the  turbulence  pattern 
helps  to  maximize  heat  transfer.  This  is  the  idea  one  gets  from  perusal,  for 
example,  of  the  papers  by  Gardon  and  associates  (References  7,9),  Chamberlain, 
Metzger  and  Walz,  (References  8,30,31)  where  heat  flow  patterns  are  studied 
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at  the  stagnation  point  with  the  Reynolds  number  as  a  parometer  and  Zp/D 
as  abscissa  (Figure  6).  This  phenomenon  is  somewhat  obscured  when  what 
actually  has  been  measured  Is  not  strictly  the  local  heat  transfer,  but  happens 
to  be  only  the  average  heat  flow.  This  occurs  naturally  with  the  heat  transfer- 
sensitive  element  larger  than  the  effective  jet  diameter;  it  can  be  seen  in 
Figure  6  of  the  paper  by  Huang  (Reference  42)  and  Figure  9  of  that  by  Walz 
(Reference  30).  Into  this  category  seem  to  fall  also  observations  of  Perry 
(Reference  40). 

In  general,  the  experiments  of  each  particular  investigator  have  been 
carried  out  over  only  a  relatively  limited  range  of  parameters,  thus  making  a 
general  comparison  difficult.  In  addition,  some  variables  apparently  could 
not  alway>  be  kept  completely  under  control,  which  led  to  what  appears  to  be 
contradictory  result*,  it  times  (Huang,  Reference  42,  as  pointed  out  by  Walz, 

Reference  30). 

Four  more  general  methods  for  measuring  of  heat  transfer  have  been  usea 
so  far:  heat  transfer  gauges  (Gardon  and  associates.  References  7,9, Cadek, 

Reference  33),  temperature  change  of  a  piece  of  metal  of  high  conductivity  in 
the  nonsteady  state  (Metzger,  Huang,  Schauer  and  Eustis  and  Walz,  References 
31,42,14,30),  metal  rod  calorimeters  (Chamberlain,  Tomich,  Hrycak  and  Comfort  et  al., 
References  8,35,39,44)  and  actual  steam  calorimeters  (Kezlos,  Reference  2). 
Additional  measuring  techniques  are  possible  with  mass  transfer  analogy  (cf. 

Koopman  and  Sparrow,  and  Schrader,  References  54,70). 

While  most  investigators  considered  heat  transfer  from  impinging  jets  to 
flat  surfaces,  Walz  (Reference  30)  concerned  himself  also  with  effects  of 
axially  symmetric  jets  Impinging  on  spheres  externally,  and  two-dimensional  jets 
impinging  externally  on  cylinders,  with  the  long  axis  of  the  jet  impinging 
parallel  to  the  cylinder  axis. 

Some  Investigators  considered  jet  Impingement  effects  on  heat  transfer 
for  angles  other  than  normal  (Perry,  Schauer  and  Eustis,  References  40  and 
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14),  and  several  Investigators  pursued  the  effect  of  cross-flow  on  heat  transfer 
from  the  Impinging  jets,  extremely  Important  for  many  application'  Heat 
transfer  resulting  from  a  single  jet  In  cross-flow  has  been  investigated  by 
Bouchez  and  Goldstein  (Reference  71). 

Wall  jet  heat-transfer  predictions,  all  being  of  semi -empirical  kind,  have 
been  found  generally  to  agree  with  experimental  results.  The  theoretical 
approaches  employ  the  consideration  of  two  flow  layers:  an  inner  layer  along 
the  wall  and  an  outer  layer  assumed  to  behave  as  a  free  jet.  Both  Integral 
and  finite  difference  solutions  have  been  obtained  by  various  Investigators 
(cf.  References  14,33,39,73). 

5.2  Arrays  of  Jets 

5.2.1  Arrays  of  round  jets  represent  a  configuration  which  is  perhaps  the  most 
applicable  for  cooling  of  large  areas.  Care  must  be  taken  to  properly  space 
the  holes  so  that  jet  Interference  is  minimized,  with  spaclngs  of  from  4 
to  6  holes  diameters  (about  5  to  2  percent  of  free-flow  area).  Results  of  tests 
reported  In  Reference  7  showed  that  the  best  correlation  could  be  obtained 
by  use  of  the  velocity  of  arrival,  and  the  center-to-center  hole  spacing  as 
the  characteristic  values,  as  a  single  jet  yielded  better- heat-transfer  results 
than  one  jet  In  an  array  of  jets.  Relatively  good  agreement  was  found  by 
comparing  the  results  of  several  different  investigations  of  arrays 
of  round  jets.  The  effect  of  crossflow  on  heat  transfer  may  be  accounted 

for  by  means  of  an  experimentally  obtained  correction  parameter;  generally, 
a  decrease  in  heat  transfer  resulted  from  Increasing  the  crossflow  to  the  jet. 
Another  way  the  crossflow  effects  could  be  accounted  for  In  determining  the 
average  Stanton  number,  by  including  In  the  correlation  a  power  of  the  ratio 
of  the  crossflow  to  the  jet  flow. 
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5.2.2  Arrays  of  Slot  Jets  -  Data  for  two  distinct  arrays  of  slot  jets  showed 
that  the  identity  of  each  jet  was  preserved  at  small  slot-to-plate  spacings, 
but  jet  Interaction  before  impingement  occurred  at  larger  spacings.  .he  average 
data  obtained  for  multiple  slots  have  been  successfully  correlated.  It  appears 
that  slot  jets  may  be  more  Immune  than  round  jets  to  the  effects  of  crossflow. 

When  a  crossflow,  equal  to  as  much  as  60  percent  of  the  air  flow,  was  <jper- 
imposed  on  an  array  of  slot  jets,  no  reduction  in  heat  transfer  was  found. 

5.3  Recommendations 

It  appears  that  a  great  amount  of  generally  reliable  information  exists 
already  on  the  use  of  the  individual  jets  (of  the  round  and  slot  type)  for  cooling, 
or  heating  of  small  areas  with  high  Intensity  heat  fluxes  (spot  cooling  or  heat¬ 
ing).  For  cooling  of  larger  flat  areas,  arrays  of  jets  appear  to  be  appropriate 
and  their  performance  under  laho-cory  conditions  can  be  described  with  a  reason¬ 
able  degree  of  ?c..urac'\  T‘.a  point  Is  made  here,  however,  that  in  actual  appli¬ 
cations,  where  arraying  of  a  large  number  of  jets  becomes  mandatory,  the  effects 
of  crossflow  may  substantially  reduce  the  effect  of  impingement  cooling.  Also 
proper  evaluation  of  the  effect  of  spent  air  becomes  here  a  problem  of  utmost 
importance. 

While  there  is  generally  satisfactory  Information  on  impingement  cooling/ 
heating  of  flat  and  concave  surfaces,  virtually  no  information  exists  on  cooling 
or  hrating  of  large  convex  surfaces.  Experimental  work  should  be  recommended, 
therefore,  In  this  area,  as  well  as  additional  theoretical  studies.  It  is  believed 
that.  In  the  near  future,  there  will  be  a  substantial  number  of  new  contributions 
In  the  area  of  Impinging  jet  heat  transfer,  of  tneoretical  and  practical  signifi¬ 
cance.  New  technological  development  will  be  generated  by  proper  organization 
and  channelling  of  research  related  to  impinging  jets  In  general,  and  to  heat 
transfer  from  impinging  jets  in  particular. 
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Potential  Flow  Solutions  at  Stagnation 
Experimental  Results. 
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FIGURE  5  Comparison  of  Poplel's  Results  with  Those  of  Other  Investigators. 


-  55  - 


FIGURE  7  Oyban  and  Mazur's  Stagnation  Point  Meat  Transfer  Measurements 
Carried  Out  on  Turbullzed  Jets. 
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FIGURE  8  Distribution  of  Dimensionless  Velocity  Gradient  at  Stagnation 
Point;  Comparison  of  Results  of  Various  Investigators. 
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FIGURE  9  Comparison  of  Heat  Transfer  Results  for  Rows  of  Round  Jets  and 
Arrays  of  Round  Jets: 

(a)  Kercher  and  Tabakoff,  (b)  Hrycak,  (c)  Gardon  and  Cobonpue 


FIGURE  10  Comparison  of  Average  Heat  Transfer  Results  for  Rows  and  Arrays 

of  Round  Jets:  -  Results  adapted  from  Reference  60  (small 

diameters  throughout). 
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FIGURE  11  Additional  Comparison  of  10  by  5  In.  Semi  cylindrical  Plate  Results 
with  Results  of  Other  Investigators. 
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FIGURE  72  Local  Distribution  of  Heat  Transfer  Coefficients,  10  by  5  in. 

Semi  cylindrical  Plate  (Reference  39);  Comparison  with  Results  of 
Other  Investigators. 


